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Abstract. In this paper we present the first results of deep star 
counts carried out within the Calar Alto Deep Imaging Survey, 



CADIS (Meisenheimer et al. 1998). Although CADIS was de 



signed as an extragalactic survey, it also attempts to identify the 
stars in the fields in order to avoid confusion with quasars and 
compact galaxies. 

We have identified a sample of about 300 faint stars (15.5 Si 
R ^ 23), which are well suited to study the structure of 
the Galaxy. The stars lie in two fields with central coordi- 
nates a 20 oo = 16 h 24 m 32. s 3, <5 20 oo = 55°44'32" (Galactic 
coordinates: I = 85°, b = 45°) and a 20 oo = 9 h 13 m 47. s 5, 
£2000 = 46°14'20" (I = 175°, b = 45°) (hereafter 16h and 9h 
field, respectively. The stars have been separated from galaxies 
by a classification scheme based on photometric spectra and 
morphological criteria. Distances were derived by photometric 
parallaxes. We are able to find stars up to distances of 25 kpc 
above the Galactic plane. The vertical density distribution of 
the stars shows the contribution of the thin disk, the stellar halo 
and the "thick disk" of the Galaxy. We give quantitative de- 
scriptions of the components in terms of exponential disks and 
a de Vaucouleurs spheroid. For the disk stars we derive the lu- 
minosity function. It is equal within the errors to the local lu- 
minosity function and continues to rise out to at least My = 13. 
Implications for the mass function are briefly discussed. 
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1. Introduction 

The stellar structure of the Galaxy has been studied by many 
groups using a variety of methods, such as metallicity-age 
determinations or studies of the kinematics of nearby halo and 
disk stars (for a review see Norris 1998). Other discussions 
can be found in Fuchs & JahreiB (1998), or Fuhrmann (1998). 
The method predominantly used to study Galactic strucure is 
the method of starcounts. This provides a measurement of the 
density distribution of the stellar component of the Galaxy 
dGilmore & Reid 1983| ; fceid et al. 199(1 fceid et al. 1997| ; 



pould et al. "l998| ). 

In the so called standard model (Bahcall & Soneira 1980) 
the vertical structure of the disk follows an exponential law 
(p ~ e" z//l2 ) with scaleheight h z . As Gilmore & Reid (1983) 
showed, the data can be fitted much better by a superposition 
of two exponentials. Whether or not this deviation from the 
single exponential is due to a distinct population of stars is not 
clear, although this is su ggested by the d ifferent kinematics 
and lower metallicities ( [Freeman 1992 ). We will call the 
deviation "thick disk", regardless of its unknown physical 
origin. 

Although the existence of the "thick disk" component seems 
well established now ( jVorris 1998| ), the proper values for 
the scaleheights and lengths are not exactly determined yet, 
since in most studies the distance of observable stars has been 
limited to a few kpc above the Galactic plane. 
Another topic of debate has been whether the stellar luminosity 
function (hereafter SLF) of main sequence stars in the disk 
declines beyond My = 12 mag . Most photometric studies find 
a down-turn of the slope o f the SLF at M v = 12 mag (Stobie 
et al. 1989; [Kroupa 1995|; |Gould et al. 1998|), whereas local 



observations of stars in the solar neighbourhood indicate a flat 
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continuation to fainter magnitudes ( Wielen et al. 1983 ; JahreiB 
& Wielen 1997a). Therefore, it is worthwhile to derive the 
SLF from the CADIS starcounts. 

The paper is structured as follows: the Calar Alto Deep 
Imaging Survey is outlined in Sect. 2. Sect. 3 describes the 
preparation of the stellar data, Sect. 4 deals with the density 
distribution of the stars, Sect. 5 with the SLF and its implica- 
tions for the mass function. Sect. 6 gives a brief summary and 
an outlook on future prospects. 



2. The Calar Alto Deep Imaging Survey 

The Calar Alto Deep Imaging Survey combines an emission 
line survey carried out with an imaging Fabry-Perot interfer- 
ometer with a deep multicolour survey using three broad-band 
optical to NIR filters and up to eighteen medium-band filters 
when fully completed. The combination of different observing 
strategies facilitates not only the detection of emission line 
objects but also to derive photometric spectra of all objects 
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Table 1. Exposure times for the three broadband filters used 
for the analysis of the stellar component in in the two CADIS 
fields. Numbers shown in italics indicate, that the exposures 
were actually obtained at the 3.5 m telescope and scaled to the 
2.2 m telescope. 



Filter t. 



exp,16h 



exp,9h 



B c 
Rc 

ZB15 



6200 10000 
5300 2800 
30700 1750 



in the fields without performing time consuming slit spec- 
troscopy. Details of the survey and its calibration will be given 
in Meisenheimer et al. (in preparation). 

All observations were performed on Calar Alto, Spain, in the 
optical wavelength region with the focal reducers CAFOS 
(Calar Alto Faint Object Spectrograph) at the 2.2 m telescope 
and MOSCA (Multi Object Spectrograph for Calar Alto) at 
the 3.5 m telescope, and with the Omega Prime camera for the 
NIR observations. 

As a byproduct of the survey we obtain a lot of multi-color 
data about faint stars in the Galaxy. Although for the ob- 
ject classification (see below) exposures in two broadband 
filters and seven medium-band filters (in the case of the 
9h field eight medium-band filters) are used, the present 
analysis of the stellar component of CADIS is based only 
on exposures in three filters, Rc (central wavelength/width 
A C /AA = 649nm/170nm), B c (A c = 461 nm/lOOnm) and 
-^815 (A c = 815nm/32nm). Exposure times converted to the 
2.2 m telescope are given in Table |lj 

The nine CADIS fields measure w 1/30 0° each and are 
located at high Galactic latitude to avoid dust absorption and 
reddening. In all fields the total flux on the IRAS 100 [im maps 
is less than 2MJy/sr which corresponds to Eb-v < 0.07, 
so we do not have to apply any color corrections. A second 
selection criterium for the fields was that there should be no 
star brighter than w I6 mag in the CADIS R band. In fact the 
brightest star in the two fields under consideration has an R 
magnitude of 15.42 mQ f . 



2.1. Object detection and classification 

Objects are identified on each of the deep images (super- 
position of 5 to 15 individual exposures) using the Source 
Extractor software SExtractor (Bertin & Arnouts 1996), and 



the resulting lists merged into a master catalogue. 
Photometry is done using the program Evaluate, which has 
been developed by Meisenheimer & Roser (1986). Variations 
in seeing in between individual exposures are taken into 
account, in order to get accurate colors. For photometric 
calibration we use a system of "tertiary" standard stars in the 
CAD IS fields, which are calib rated with secondary standard 
stars ( Oke 1990 ; Walsh 1995 ) in photometric nights. 
From the locus of the stars in the 8 -dimensional color space 



we conclude that the relative calibration between each pair of 
wavebands is better than 3 % for all objects with R = 22. 
Since one of the major goals of the survey is the classification 
of every object found in all CADIS fields (w 80 000 to 100 000 
in total), a classification scheme was developed which is based 
on template spectral energy distributions (see Wolf 1998). The 
observed colors of every object are compared with a color 
library of known objects, whose colors are obtained from 
synthetic photometry performed on our CADIS filterset. The 
input library for stellar spectra was the Gunn & Stryker (1983) 
catalogue. For each object the probability to belong to a certain 
object class (stars - quasars - galaxies) is computed. 
Objects classified as stars have stellar colors with a likelihood 
of more than 75%, and images the profile of which does not 
deviate significantely from that of well defined stars. 
Details about the performance and reliability of the classifica- 
tion are given in Wolf et al. 1999, and Wolf (in preparation). 
With the current filter set and exposure times the classification 
is reliable down to a limit of R ~ 23 maff . This was checked 
by spectroscopic follow-up observations of 245 arbitrarily 
chosen objects (55 stars, 153 galaxies and 20 quasars) with 
R < 23 ma9 . One galaxy has been classified as a star by it's 
colors, two quasars as galaxies and two galaxies as quasars. 
The star counts are not significantly affected by the misclas- 
sifications down to R = 23 maff . Thus we restrict our present 
analysis to stars with R < 23 mag . 



3. Preparation of the stellar data 

3.1. Photometric parallaxes 

To derive distances of the stars from the distance modulus 
mji — Mb,, it is essential to know the absolute magnitudes of 
the stars, Mr. 

In principle absolute magnitudes of main sequence stars can 
be obtained from a color-magnitude diagramm. This main se- 
quence approximation is valid for all stars in our sample since 
we can be sure that it is free from contamination of any non- 
main sequence stars, as the faint magnitude intervall we ob- 
serve (16 < R < 23) does not allow the detection of a giant 
star. 

We took a mean Mvj versus (B — V)j relation from Lang 
(1992). A complication arises, because there is no V filter 
included in our filter set. Thus we have to convert My and 
(B — V)j into our filter system, in order to derive the absolute 
magnitudes from mean main sequence fit in Mr versus (b — r), 
where the CADIS color is defined by: 



(6-r)=2.5Iog F7 



(1) 



Here F£ is the flux outside the atmosphere in units of Photons 
m~ 2 s _1 nrn -1 in the CADIS filter k. 

In order to convert (b — r) into the Johnson (B — V)j we per- 
formed synthetic photometry on the Gunn-Stryker stars (Wolf, 
priv. comm.). 
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b — r can be calibrated to the Johnson-Cousins system (the 
CADIS Re is very close to the Cousins R) by using Vega as a 
zero point: 

^ CS a(A = 440nm) 



KLa(A = 648nm ) 



(B-R) c = (6 -r) + 2.51o; 

= (6-r) + 0.725. 

The absolute magnitude Mr c is then given by 

M Ro = M Vj - (Vj - R c ) , with 
(Vj - R c ) = (B c - Rc) ~(B- V)j , 



(2) 



(3) 



where we assume Bj = Be (see Huang et al., in preparation). 
With the above conversions, the main sequence (Mr vs (b—r)) 
can be approximated by a fourth order polynomial in the range 
-1 < (b-r) < 1.8: 



M R = c + ci(b-r)+ c 2 (b - rf 
+ c 3 (b-r) 3 +c 4 (b-r) 4 , 

the parameters of which are: 



(4) 



c = 4.01236 ci = 4.12575 c 2 = -1.89076 
c 3 = 0.762053 c 4 = 0.341384, 

as shown in Fig. |jj 

One further complication arises due to the fact that the mean 



sequence relation, but shifted towards fainter magnitudes by 
Mr = 0.75 ma9 (see Fig. [!]). This value is the mean devia- 
tion from the mean main sequence defined by the CNS 4 stars 
( Jahreifi & Wielen 1997b[ ) of a subsample of 10 halo stars, 
for which absolute R magnitudes were available (Jahreifi, priv. 
com.).| 

The spread in a two color diagram (b — r) versus (r — i) (that 
is the CADIS color between Rc and Jgis analog to Eq. (1)) 
becomes significant at (b — r) ~ 1.0, see Fig. ^. The maxi- 
mal photometric error for the very faint stars is 0.15 ma9 . Here 
metallicity effects will distort the relation between the mea- 
sured (b — r) colors and the spectral type (temperature) and 
thus lead to wrong absolute magnitudes, so we have to correct 
for metallicity in order to avoid errors in the photometric par- 
allaxes. 

The Rc filter is strongly affected by metallicity effects like 




-1 o 1 

(b-r) 



Fig. 1. The mean main sequence from Lang (1992), converted 
to CADIS color b — r and Rc magnitude; the solid line is the 
fourth order polynomial (Eq. (Q)). For the blue (b — r < 0.7) 
metal-poor stars the mean main sequence is shifted towards 
fainter magnitudes. 

main sequence relation is valid strictly only for stars with so- 
lar metallicities, whereas our sample may contain stars spread 
over a wide range of different metalicities. For the blue stars 
(b — r < 0.7), which are almost all halo stars (see Fig. Q), and 
therefore supposed to be metal-poor, we use the same main 



0.5 1 1.5 2 

(b-r) 

Fig. 2. The b—r colors are projected along the "isophotes" 
onto the mean main sequence track (solid line). For some of 
the very red stars (b — r > 1.8), which are subject to the largest 
corrections, photometric errors are plotted. 

absorption bands of Ti02 and VO molecules in the stars' atmo- 
sphere, whereas the Be and the medium-band filter I^is (the 
wavelength of which was chosen in order to avoid absorption 
bands in cool stars) are not. So in a first approximation we can 
assume the "isophotes" of varying metallicity in a (b — r) ver- 
sus (r — i) two color diagram to be straight lines with a slope 
of —1, along of which we project the measured colors with 
(b — r) > 1.0 onto the mean main sequence track which in the 
interval 1.0 < (b — r) < 1.8 is defined by 

(r-i) = 0.39 (b - r)i mi - 0.36 (b - rf COIt 

+ 0.09 (b - r) + 0.06 . (5) 

This projection implies that stars with (b — r) COII ^1.8 cannot 
exist in Fig. ^, which shows the spatial distribution of metal- 

1 Note that this artificial separation may well lead to wrong absolute 
magnitudes for individual stars (e.g. a disk star with r < 1 kpc but 
b — r < 0.7), but should be correct on average. 
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Fig. 3. Spatial distribution of metallicity corrected (b — r) corr . The solid line represents the distance dependent upper color limit 
at R = 23 mag , the dotted line is the lower limit due to the selection criterium of the fields (no star brighter than R as 15.5). 
Since the metal-poor halo stars are intrinsically fainter, the color limits are shifted accordingly. The horizontal line denotes the 
cut between halo and disk stars, the dashed-dotted line indicates the cutoff at (b — r) corr = 1.8 due to the metallicity correction. 
The different symbols refer to the different fields: squares - 16h field, triangles - 9h field. 



licity corrected (b — r) corr colors (the limit is indicated by the 
dashed-dotted line). Both the upper and lower magnitude lim- 
its lead to selection effects which have to be taken into account. 
As one can see in Fig. || there is an bimodality of the observed 
color distribution. The accumulation of red stars in the upper 
left consists mainly of disk stars, and is separated by a void 
from the blue stars, which predominantly belong to the halo. 
Thus a crude disk-halo separation can be drawn by a color cut 
- we take stars with (b — r) corr < 0.7 to be halo, stars with 
(b — r) COII > 0.7 to be disk stars. In the following we will 
make use of this color cut to derive the distribution of the disk 
stars separately. In a second step, the distribution as a whole 
will be analysed. In Fig. |]the cut is denoted by a dotted hori- 
zontal line. 

In the two fields we have analysed so far we find 95 halo and 
178 disk stars, that is a factor of two more disk stars than pre- 
dicted by the standard model (Bahcall & Soneira 198C). This 
surplus of disk stars was already noted by Reid & Majewski 
(1993). Fig. ^ shows the distribution of the (b — r) corr in the 
two fields under consideration. 

4. Density distribution of the stars 

4.1. Completeness correction 

As expected, the detection limit of absolute magnitudes (col- 
ors) is distance dependent (see Fig. fy. 



16H field 
9H field 



Fig. 4. Distribution of (b — r) corr in the two fields, this bimodal 
distribution can be used to separate halo and disk stars by a 
color cut. 



We use the two-dimensional distribution of stars in the (b — 
?")corr vs l°g r diagram (Fig. |3|) to correct for this incomplete- 
ness in the following way: First, we divide the distance in log- 
arithmic bins of 0.2 as indicated in Fig. [}] and count the stars 
up to the the upper color limit (b — r) i (this is the distance 
dependent color (luminosity) limit, up to which stars can be de- 
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tected. The metal-poor halo stars are intrinsically fainter (see 
paragraph 3)), thus the limits are shifted accordingly). 
The nearest bins (—0.8 < logr < —0.2) are assumed to be 
complete. For the incomplete bins we multiply iteratively with 
a factor given by the ratio of complete to incomplete number 
counts in the previous bin, where the limit for the uncorrected 
counts is defined by the bin currently under examination (j): 



/v corr 

3 



N 



na + &>. 



1=1 



N 



(6) 



where Nj is the number of stars in bin j, N is the number of 
stars in the previous bin (J — 1), up to the limit given by the bin 
j, and N is the number of stars from that limit up to the limit 
given by bin j — 1, see also appendix. 

With the poissonian errors cat = VW, ct n > — %/iV 7 , and 
\l N" the error of the corrected number counts be- 



comes: 



2 




'N'. 



(7) 



For a detailed deduction see appendix. 

The completeness correction is done for each field separately. 
With the corrected number counts the density in the logarith- 



mic spaced volume bins (Vj — |w(r| +1 
calculated according to 



lycorr 
3 

V 



Vj 



■ )) can than be 



(8) 



For every logarithmic distance bin we use the mean height 
z above the Galactic plane < Zj >— s'mb- < r >, where 
< logr >= logr-, + (logrj+i - logrO/2 = logr, + 0.1 =^< 
r >= 1.259tv. 



4.2. Vertical density distribution 

We first study the density distribution of the disk stars by tak- 
ing only stars in the corresponding color interval into account. 
Although the color-cut at (6— r) corr = 0.7 is a rather crude sep- 
aration between disk and halo, we gain a clearer insight into the 
disk distribution since the contamination by halo stars is sup- 
pressed considerably. 

As the nearest stars in our fields have still distances of about 
200 pc the normalization at z = has to be established by 
other means: We take stars from the CNS4 (J ahreifi & Wielen 
1997b), which are located in a sphere with radius 20 pc around 
the sun. 



The stars in our normalization sample are selected from the 
CNS4 by their absolute visual magnitudes, according to the 
distribution of absolute magnitudes of the CADIS disk stars 
(6.5 < M v < 14.5). 

Fig. g shows the resulting density distribution of the disk stars 
in the two CADIS fields. The solid line represents a fit with a 
superposition of two exponentials, the dotted line is the fit for 
the thin disk component (the first seven data points). Obviously 
a single exponential is not a good description. 
It was suggested to fit the thin disk with a secans hyperbolicus 
- the exponential is unphysical in that sense, that it is not con- 
tinuously differentiable at z = 0. A squared secans hyperboli- 
cus (which represents a self-gravitating isothermal disk) can be 
proved not to fit the data very well - indicating that the stellar 
disk is in no way isothermal (the velocity dispersion depends 
on the spectral type). 

The fits for the three functions under consideration: 

Pexp{z) = Tii exp(-z//ii) + n 2 cxp(-z//i 2 ) (9) 
Psechiz) = n 3 scch(-z/z ) + n 4 cxp(-z//i 3 ) (10) 
Psechz(z) = n 5 sech 2 (-z/z ) + n 6 cxp(-z//i 4 ) . (11) 

are shown in Fig. 6, the corresponding parameters are given in 
Table g. 

As the secans hyperbolicus is a sum of two exponentials 

2 



sech(z/z ) 



cxp(z/z ) + cxp(-z/z ) 



(12) 



zq is not really a scaleheight, but has to be compared to hi by 
multiplying it with arcsech^ w 1.65745: h' x — z ■ 1.65745 
(cf. Tab. g). 

The errors of the scaleheights and corresponding parameters 
Table 2. Fit parameter for the three fit functions 



Parameter 


16h field 


9h field 


m 


5.06-10 7 


5.84 -10 7 


"2 


4.14 -10 6 


6.23 -10 5 


hi 


(280 ± 14)pc 


(267 ± 9)pc 


h 2 


(1267 ± 74)pc 


(1296 ± 116)pc 


713 


4.54 -10 Y 


5.19 -10 Y 


"4 


9.86 -10 6 


1.02 -10 6 


«0 


(172 ± 12)pc 


(216 ± 8)pc 


hi 


(286 ± 20)pc 


(358 ± 13)pc 


h 3 


(1478 ± 30)pc 


(1132 ± 84)pc 


n 5 


4.62 -10 Y 


5.4 -10' 


«ii 


2.78 -10 6 


1.37 -10 6 


zo 


(352 ± 26)pc 


(360 ± 13)pc 


fl4 


(1158 ± 112)pc 


(1047 ±71)pc 



zq and zq are estimated by changing their values until \ 2 
increases by 1. We find that at z — the contribution of the 
thick disk component to the entire disk is 2 ± 4%. 
Our values for the scaleheights lie within the range given by 
different authors in the literature. The proposed values for the 
exponential scaleheights of the thin disk range between 200 
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/kpc 



r > 
data 



Fig. 5. Vertical density distribution of the disk stars (b 
0.7) in the two fields, the normalisation (star) is given by 
from the CNS4. The distribution shows clearly the contribu 
tion of the thick disk component - the solid line is a fit v 
superposition of two exponentials, the dotted line is the 
exponential fit for the thin disk component only. 



sin 



pc flKentetal. 1991 [ ) and 325 pc ( fcahcall & Soneira 1 
The parameters £q and /14 we found in the 16 h field are 



within the errors to the values found by Gould et al (1997 ) 



4.3. Density distribution in the halo 



Fig. [7| shows the density distribution of all stars in the 
CADIS fields. The dotted line is the secans hyperbolicus 
ponential fit for the thin and thick disk components, the d, 
line is a deVaucouleurs law. Although the space density 
sponding to the r 1 ' 4 law in projection has no simple am 



two 
i jt- ex- 
: shed 
c orre- 
lytic 

form ( jYoung 1976| ) it is possible to form expansions abot t the 

1 for 
dBah- 



origin and the point at infinity. The latter expansion is 
a large range of distances. An analytic approximation is 



val:d 



ith a 
gle 



)80) 



< :qual 




/kpc 



, A secans hyperbolicus (dotted line) or a squared secans 
hyperbalicus (dashed line) can be shown to fit the thin disk 
well like an exponential (solid line). The normalisation 
1 by data from the CNS4. 



call H86): 



1.21 i 



where 
I are C 
from 
x = 
dot 
tion 
The 
with th 
The 
ities 
selecte 
((6-r 
Fig _ 
to z = 
couleu 
disk fit 
« 1.5 



ar d 



contrib ution 
disk arjd 
If we 
ter and 
CNS4, 



exp 



z,bj) = po- 



10. 093 



(4>) 



1/4 



10.093 



7/8 



exp 



-10.093 



(*) 



1/4 



10.093 



6/8 



(*)' 

0.08669/ {r/r q ) 1/4 } ,r> omR Q 



,R< 0.037? 







(13) 

-4t - 2i? -^-cosZ) 1/2 , b and 

tan* 2 b u tan b > ' 

alactic latidude and longitude, the distance of the sun 
Galactic center R Q = 8kpc, R = (x 2 + z 2 ) 1/2 , 
I + d 2 cos 2 b — 2R Q d cosb cos Z) 1 / 2 , z = dsinb. The 
dashed is the sum of both halo and disk distribution func- 



R = (R 2 + z 



t le 



dejnsity of halo stars extrapolated to z = is compared 
data from the CNS4 with an estimated error of w 10%. 
staYs are discriminated against disk stars by their metallic- 
kinematics (Fuchs & JahreiB 1998). They are further 
1 according to the color cut of the CADIS halo stars 
< 0.7). 

ihows the distribution of all stars in the range z = 
10 kpc, the fit for a single exponential disk, the deVau- 
s law, the sum of both, and the sum of the thin plus thick 
and the deVaucouleurs law. It is obvious that between 
kpc to 5.0 kpc the thick disk provides the predominant 
to the overall distribution. Thus, the sum of thin 
deVaucouleurs halo does not fit the distribution, 
issume the sun to be 8 kpc away from the Galactic cen- 
the normalisation po to be fixed by the data from the 
all parameters are completely determined. As can be 
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1 



\ 



/kpc 




z /kpc 

Fig. 7. Density distribution perpendicular to the Galactic Plane, 
stars are data taken from the CNS4; the dotted line is the secans 
hyperbolicus plus exponential fit for the disk, the dashed line is 
a deVaucouleurs law and the dot-dashed the sum of both. 

seen in Fig. ^, the corresponding plots fit the distribution of the 
halo stars up to a distance of about 22 kpc above the Galactic 
plane. 

The de Vaucouleurs law is an empirical description of the 
density distribution s of stars in ellipticals, bulges of spirals 
and galactic halos ( ]de Vaucouleurs & Pence 1978 ) but it is 
equally justified to fit other distribution functions to the data: a 
smoothed power law fits the halo stars equally well (Fuchs & 
JahreiB 1998; |Gould et al. 1998| ), see e.g. Fig. § When apply- 
ing the power-law fit, we add a further free parameter which 
describes the flattening of the halo: 



PH{z,b,l) 



a/2 



J c ~r '© + ta^b 



(14) 



with an arbitrarily chosen core radius of r c = 1 kpc. 

We find an axial ratio c/a of (0.63 ± 0.07), with an exponent 

a = (2.99 ±0.13). 



z /kpc 

Fig. 8. Distribution of all stars in the range << z < 10 kpc, 
the fit for a single exponential disk (dotted line) the deVau- 
couleurs law (short-dashed), the sum of both (dot-dashed), and 
the sum of the thin plus thick disk fit and the deVaucouleurs 
law (long-dashed). 



9h field 



l/r<" 

de Vaucouleurs 



z /kpc 

Fig. 9. An inverse power law (broken line) can also be fitted to 
the data - in comparison with the a de Vaucouleurs law (solid 
line), shown here for the CADIS 9h field. 



5. The stellar luminosity function 

The knowledge of the distribution function enables us to calcu- 
late the SLF for the thin disk stars. We selected the stars with 
distances less than 1.5 kpc. Beyond this point the contribution 
by thick disk and halo stars becomes dominant. 
Although CADIS does not include a filter which is close to 
the Johnson-V, we calculate absolute visual magnitudes from 
the R magnitudes to make comparison with literature eas- 
ier. In the relevant magnitude intervall (5 < My < 14, i.e 
— 1 < (6 — r) < 1.8) there holds a linear relation (see Sect. 
O): 



M Vj = 1.058 • M Rc 



(15) 
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We calculate an effective volume for every luminosity 
integrating the distribution function along the line of 
where the integration limits are given by the minimum 
1.5 kpc and the distance modulus derived for upper 
limiting apparent magnitude: 



V* K = w I v(r,b)r*dr, 



a id 



where 

d _ in 0.2(16 mas -A/ ii )-2.0 

i? max = min(1.5kpc,10 02 ( 23mas - M «)- 20 ). 
The distribution function 

v{r,b) = exp(— r sin6//ii) 

is normalised to unity at z = 0. 
The weighted mean of the SLFs of the two CADIS 
shown in Fig. |l^, in comparison with the SLF of the 



side a distance of 20 pc ( Jahreifi & Wielen 1997a ) 
based on HIPPARCOS parallaxes. As can be seen f 



(17) 

fields is 
stars in- 
which is 
om Fig. 
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Local Luminosity Function 
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Fig. 10. The m ean SLF of the two CADIS fields (solid 
the local SLF (ttahreiB & Wielen 1997ah. 



[l Oj , the CADIS SLF is equal within the errors to the kjcal SLF. 
Since the errors are dominated by Poissonian statistic >, the er- 
ror bars are small at the faint end, and larger at the br ight end, 
complementary to the local SLF. 
Thus the weighted mean of CADIS and local SLF 
bines the large number of bright stars in the local 
our superior sampling of faint stars (shown in Fig 
regarded as the most accurate determination of the 
combined CADIS/local SLF keeps rising with constan 
its limit at My = 13. This is in pronounced contrast 
determinations of the SLF based on faint star counts 
al. 1989; Kroupa 1995). The discrepancy is demon 
Fig 



san pi 



i ( Stobie et 
s trated in 
:he most 

recent photometric SLF based on HST observations ( Sould et 



1 1 where we compare the combined SLF with 



bin by 
sight, r, 
between 
lower 



(16) 



line) and 



wh|ch com- 
e with 
can be 
$LF. The 
slope to 
tolprevious 



flint 



al. 1998|) 
At this 

exact lodation 
errors wl lich 



Fig. 11. 

parison 
vations 



' The weighted mean of CADIS and local SLF in com- 
\ /ith a photometric SLF which is based on HST obser- 
( Zrould et al. 1998|). 



5.7. Imp 

The stel 
from the 
stars (r 
relation 
Henry e 
masses f 
in equally 
counts 
(see Eq 
the SMF 



end both incompleteness and uncertainties of the 
of the main sequence may introduce systematic 
are hard to quantify. 



Weighted mean of CADIS and local LF 
HST 
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ications for the mass function 



ar mass function (SMF) can in principle be inferred 
SLF by converting the My magnitudes of the disk 
1.5 kpc) into masses; this requires a mass-luminosity 
for which we adopted the analytical fit taken from 
al. (1993, 1999). Their relation is valid for stellar 
-om 0.08 M /M Q to 2.0 M/M & . If we count the stars 
spaced mass intervals of 0.1 and divide the number 
the jth interval by the maximum effective volume 
[To])), calculated for the corresponding luminosity, 
can be represented by a power law: 



o) <* 



M 



(18) 



From thf data of the disk stars in the two CADIS fields we 
derived 
is equal 
al. (19931) 



= -1.28±g;|| for [0.2 < M/M Q < 1.1]. Its slope 
within the errors) to the value proposed by Henry et 
(o = -l). 



6. Sumnpary and future prospects 

Although the Calar Alto Deep Imaging Survey is designed 
as an e: :tragalactic survey, we obtain a substantial set of 
multicoL >r-data about faint stars in the Galaxy. With the cur- 
rent filte • set and exposure times the classification is reliable 
down to \Mfi = 23. 
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From the w 300 stars we identified in two CADIS fields cover- 
ing a total area of 1/15 □", we deduced the density distribution 
of the stars up to a distance of about 20 kpc above the Galac- 
tic plane. The density distribution shows unambiguously the 
contribution of a thick disk component of the Galaxy with a 
scaleheight of sw 1.3 kpc. 

There has been discussion whether this "thick disk" is intro- 
duced artificially by the assumption that all stars are on the 



main sequence (Bahcall & Soneira 1984). Our present sample 



of very faint field stars does rule out this ambiguity: if there 
were any giant stars included in our data, they would have dis- 
tances of about 250 kpc, thus their contribution to our sample 
can be neglected. 

The density distribution in the halo, which essentially is cor- 
roborated by comparing the local density of CNS4 stars with 
the density between 5.5 kpc and 15 kpc above the Galactic 
plane from CADIS, is perfectly fit by a de Vaucouleurs' law. 
Based on our present, limited data set of 72 stars beyond z=5.5 
kpc, it is equally well justified to fit an inverse power law with 
exponent a = 2.99 ± 0.07 to the data, with the axial ratio of 
the halo as an additional free parameter. The best fit value for 
the axial ratio is (c/a = 0.63 ± 0.07). 

Based on this knowledge of the density distribution, we deter- 
mined the stellar luminosity function (SLF) for the disk stars. 
To this end, we confine our sample to stars with distances less 
than 1.5 kpc, for at this point the contribution of the thick disk 
and the halo population becomes predominant. Our result is 
within the errors indistinguishable from the local SLF (JahreiB 
& Wielen 1997a). Thus we conclude that the weighted aver- 
age between the local SLF (containing predominantly stars at 
bright magnitudes (My < 10)) and the CADIS SLF (with su- 
perior statistics at the faint end (My > 10)) can be regarded 
as the best estimate of the SLF of the disk stars. The SLF con- 
tinues to rise at least up to My — 13, in contrast to most other 
photometric SLFs, which show a down-turn at My = 12. 
We regard this as a hint that incompleteness corrections and the 
closely connected issue of the "true" sampling volume of these 
faint star counts might need some critical revision. 
The stellar mass function (SMF) of the disk stars was derived 
by convert ing the visual mag nitudes into masses (Henry et 
al. 1993; ^enry et al. 1999| ). For a power-law SMF * oc 



{M/MqY we find a slope a = -1.28±g;||. This is (within 
the errors) consistent with the slope proposed by Henry et al. 
(1993) (a = -1). 

The main purpose of the present paper was to explore to what 
extent the faint star counts in CADIS can be used for deter- 
mining the stellar density distribution and the stellar luminos- 
ity function of the disk. However, one should keep in mind that 
our present analysis is based only on about a quarter of of the 
entire Calar Alto Deep Imaging Survey data. In the entire sam- 
ple which will contain eight fields of 1/30 □°each, we expect 
to identify a total of ^ 1200 stars with R < 23, of which 800 
are supposed to be disk stars and 400 to be halo stars. 
In addition with the complete multicolor data of the survey it 
should be possible to push the limit for a reliable star classifi- 
cation to R = 23.5 or beyond. 



As soon as better template spectra of stars with sub-solar metal- 
licities become available, the currently rough correction for 
metallicity effect can be replaced by a more accurate metal- 
licity dependent main sequence. 

With this much more accurate determination of photometric 
distances and an at least 4-fold increased statistics it should be 
possible to adress the following issues concerning the density 
and luminosity function of the stars in the Galaxy: 

- the scalelength h r of the density distribution of stars in the 
disk. This can be done if one takes the longitude dependent 
radial decrease of the density into account, e.g. the scale- 
height measured in the CADIS 16h field (I « 90°) should 
be only very slightly affected by the radial density gradient, 
whereas the 9h field is suffering from a maximal radial de- 
crease in density, so one measures an effective scaleheight 
h e fi. Likely the star counts in the 18 h and 23 h field (at 
Galactic coordinates b = 30°, I = 95°, and b = -43°, 
I = 90°, respectively, are only very slightly affected by the 
radial density gradient, so it should be possible to measure 
the "true" scaleheights hi and hi for thin and thick disk 
component with high precision. From this, the counts in 
the 9h field and the lOh field (b = 53°, I = 150°), we will 
deduce the scalelength r*j = [l/h e g — l/h\~ l . 

- the axial ratio of a flattened halo and the exponent of the 
inverse power law with higher precision. 

- a bulge - halo separation. By analysing the star counts in all 
eight fields it will be possible to separate the contribution 
of the bulge from the halo distribution. 

- the slope of the SLF with higher precision. From the entire 
data set with limiting magnitude R > 23.5 we will be able 
to deduce the SLF down to My > 14 with good statistics 
and thus decide wh ether the features like the Wielen dip 
( Wielen et al. 1983 ) at My w 7 and the apparent rise of 
the SLF beyond My — 12 are real. 



Acknowledgements. The authors would like to thank H.JahreiB who 
identified the subset of stars from the Forth Catalogue of Nearby Stars 
(CNS4) and thus provided the local normalisation for the disk compo- 
nent. 

We are greatly indebted to the anonymous referee who pointed out 
several points which had not received sufficient attention in the origi- 
nal manuscript. This led to a substantial improvement of the paper. 
We also thank the Calar Alto staff for their help and support during 
many observing runs at the observatory. 



10 



S. Phleps et al.: CADIS deep star counts 



Appendix A: Completeness correction 

Since the first distance bin is assumed to be complete, we can 
iteratively correct the number counts in the next bins. 
Nj is the number of stars in bin j, N is the number of stars in 
the previous bin (j — 1), up to the limit given by the bin j, and 
N is the number of stars from that limit up to the limit given er 



'n"" 



N 2 N ± 



and for the third bin 



by bin j — 1, as indicated in Fig. A. 1 



The division into N and N is done to avoid correlated errors. 
The first bin (j = 1) is assumed to be complete (N{). The 



NV 



N3 



Fig. A.l. Schematic sketch of Fig. |3| 



log r 




Thus the error of the completeness corrected counts in the j 
bin is 



lh 




N. 



(A.2) 



corrected number iV^ " in the second bin is then 



7V 2 corr = N 2 ■ 



Ni 



N, 



N1+N1 



N 2 \l + 



n: 



and for the third bin 



7V 3 corr = N; 




Thus this can generally be written as 



3 



N 



(A.l) 



The errors of the number counts Nj, Nj-i and iV"-_ 1 are in- 
dependent of each other and thus Poissonian: o^r 



\/Nj, 

a N ' = y 'Nj_ v a N 'i = \JN'-_ V Following Gaussian er- 
ror propagation the error for the corrected number of stars in 
the second bin is 



'N 2 



N2 
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